Systematic Effective Field Theory Investigation of Spiral Phases in Hole-Doped 
Antiferromagnets on the Honeycomb Lattice 



00 

o 
o 

(N 
D 

a^ 

(N 



F.-J. Jiang,i'0 F. Kampfer,2.| C. P. Hofmann,3.[t] and U.-J. Wiesei'4.[| 

^ Center for Research and Education in Fundamental Physics, 
Institute for Theoretical Physics, Bern University, Sidlerstrasse 5, CH-3012 Bern, Switzerland 
'^Condensed Matter Theory Group, Department of Physics, 
Massachusetts Institute of Technology (MIT), 77 Massachusetts Avenue, Cambridge, MA 02139, U.S.A. 
^ Facultad de Ciencias, Universidad de Colima, Bemal Diaz del Castillo 340, Colima CP. 28045, Mexico 
"^Institute for Theoretical Physics, ETH Ziirich, Schafmattstrasse 32, CH-8093 Ziirich, Switzerland 

Motivated by possible applications to the antiferromagnetic precursor of the high-temperature 
superconductor Naa;Co02-yH20, we use a systematic low-energy effective field theory for magnons 
and holes to study different phases of doped antiferromagnets on the honeycomb lattice. The 
effective action contains a leading single-derivative term, similar to the Shraiman-Siggia term in the 
square lattice case, which gives rise to spirals in the staggered magnetization. Depending on the 
values of the low-energy parameters, either a homogeneous phase with four or a spiral phase with 
two filled hole pockets is energetically favored. Unlike in the square lattice case, at leading order the 
effective action has an accidental continuous spatial rotation symmetry. Consequently, the spiral 
may point in any direction and is not necessarily aligned with a lattice direction. 
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I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity in the cuprates [l| , identifying the dynamical mech- 
anism behind it remains one of the biggest challenges in 
condensed matter physics. It has been suggested that 
the physics of high-tcmpcrature superconductivity can 
be described by i- J-type models. Using a variety of tech- 
niques, numerous interesting properties of doped antifer- 
romagnets have been investigate_d in great detail both 
numerically and analytically 
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For instance, as was first pointed out by Shraiman and 
Siggia a spiral phase with a helical structure in the 
staggered magnetization is a candidate ground state of 
doped antiferromagnets even at arbitrarily small doping 
Ho, 18, 21, 23, 24, 25, 26, 27, 28, 30, 31, 34, 35, 36, 37,^ 
|40|. I4I1]. Unfortunately, due to the strong electron corre- 
lations in these systems, most analytic results suffer from 
uncontrolled approximations. Similarly, numerical simu- 
lations suffer from a severe sign problem away from half- 
filling. Consequently, although numerous investigations 
have been devoted to understanding the spiral phases in 
doped antiferromagnets, some controversial results have 
been obtained. 

In analogy to chiral perturbation theory for the pions 
in QCD [12, 13 ; ^ systematic low-energy effective field 
theory for the magnons in an antiferromagnet was devcl- 
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oped in 0, i, [H, [H, [H, [H, [3. Motivated by the suc- 
cess of baryon chiral perturbation theory for pions and 
nucleons [H, [13, S, |4§|, respecting the symmetry 
constraints of the underlying t-J model and taking into 
account the location of the hole or electron pockets in 
momentum space, low-energy effective field theories for 
magnons and holes or electrons have been constructed 
for lightly doped antiferromagnets on the square lattice 
in [50, \Em, Ell • The effective theories are universally ap- 
plicable and yield results that are exact, order by order 
in a systematic low-energy expansion. Material-specific 
properties enter the effective Lagrangian in the form of 
a priori undetermined low-energy parameters, like the 
spin stiffness or the spinwave velocity c. The effec- 
tive theories for hole- and electron-doped systems were 
used to investigate the one-magnon exchange potentials 
and the resulting two-hole or two-electron bound states 
as well as possible spiral phases [Ell [H, HI, H^l • In the 
hole-doped case, the leading order magnon-hole coupling 
is described by the Shraiman-Siggia term that contains 
just a single spatial derivative. For sufficiently small ps, 
even at arbitrarily small hole density, this term stabilizes 
a zero degree spiral phase in which the spiral is oriented 
along a lattice axis. In the electron-doped case, on the 
other hand, the Shraiman-Siggia term is forbidden by 
the symmetries, and, consequently, spiral phases are not 
energetically favorable. 

In addition to the cuprates, another superconducting 
material, Naa;Co02-yH90 ISSl . has attracted a lot of at- 
tention [5^, [53, m, m, [13, l6l| . The underlying triangular 
lattice of this geometrically frustrated material leads to a 
severe sign problem and thus prevents us from studying 
it from first principles using Monte Carlo calculations. 
On the other hand, the honeycomb lattice structure of 
the dehydrated variant of Naa;Co02-yH20 at a; = 1/3 
has motivated several investigations of the antifcrromag- 
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nctism as well as the singlc-holc dispersion relation on the 
non- frustrated honeycomb lattice [G^, [sl]. In particular, 
the low-energy parameters of the effective theory for the 
t-J model, namely the staggered magnetization M.s [6^ . 
the spin-stiffness ps, the spinwave velocity c, and the ki- 
netic mass of a hole M' have been determined with high 
precision using an efficient cluster algorithm [g^l- 

Motivated by possible applications to Na,rCo02-yH2 0i 
using the same methods as for the square lattice [50, ll^l , 
we have constructed a systematic effective field theory for 
the t-J model on the honeycomb lattice. The details of 
this construction will be presented in a forthcoming pub- 
lication [gHI ■ In this work, we apply the resulting effective 
Lagrangian to investigate possible spiral phases of lightly 
hole-doped antiferromagnets on the honeycomb lattice. 
In contrast to the square lattice case, the leading terms 
of the effective Lagrangian have an accidental continuous 
rotation symmetry. This implies that possible spirals are 
not necessarily aligned with a lattice direction. Assum- 
ing that the 4-fermion couplings between holes can be 
treated pcrturbatively, the effective theory predicts that, 
depending on the values of the low-energy parameters, 
either a homogeneous phase with four or a spiral phase 
with two occupied hole pockets is energetically favored. 

The rest of this paper is organized as follows. In section 
|TT]we review the effective theory for magnons and holes in 
an antiferromagnet on the honeycomb lattice. In partic- 
ular, we list the transformation properties of magnon and 
hole fields under the symmetries of the underlying micro- 
scopic t-J model, and we discuss the accidental spatial 
rotation invariancc of the leading terms in the effective 
Lagrangian. In section Hill we consider the homogeneous 
and possible spiral phases restricting ourselves to config- 
urations that induce a homogeneous background field for 
the doped holes. In section ITVl we include the 4-fermion 
couplings using perturbation theory and investigate the 
stability ranges of the various phases. Finally, section fVl 
contains our conclusions. 



II. SYSTEMATIC LOW-ENERGY EFFECTIVE 
FIELD THEORY FOR MAGNONS AND HOLES 

In this section we briefly review the effective theory for 
magnons and holes in an antiferromagnet on the honey- 
comb lattice. In particular, we list the symmetry trans- 
formation rules for magnon and hole fields under the var- 
ious symmetries of the underlying t-J model which is es- 
sential for constructing the effective Lagrangian. The 
staggered magnetization of an antiferromagnet is de- 
scribed by a unit-vector field 

e[x) = [sin9{x)cosip{x),s\n9{x)suYLp{x),cos9{x)), (1) 

in the coset space SU{2)s/U[l)s = S"^, with x = 
(xi,a;2,i) denoting a point in (2 -I- l)-dimensional space- 
time. A key ingredient for constructing the effective field 
theory is the nonlinear realization of the global SU{2)s 
spin symmetry which is spontaneously broken down to 



its U{l)s subgroup [H^l- This construction leads to an 
Abelian "gauge" field and to two vector fields w^(a;) 

which are "charged" under U{\)s spin transformations. 
The coupling of magnons to holes is realized through a 
matrix-valued anti-Hermitean field 

Vf,{x) = ivl{x)aa, v^{x) = vl{x) T (2) 

which decomposes into an Abelian "gauge" field v^^{x) 
and two vector fields v'^{x) "charged" under the un- 
broken subgroup U{\)s- Here a arc the Pauli matrices. 
These fields have a well-defined transformation behavior 
under the symmetries which the effective theory inherits 
from the underlying microscopic t-J model 

SU{2), : v^{xy = h{x){v^{x) + d^)h{x)\ 
Di : ^'Vf^{x) = Vf^{x), 

O : = r(Ox)(iwi(O.T) + ^V2iOx) 

+ ^di + ^d2)T{Ox)\ 
^V2{x) = t{Ox){- ^ViiOx) + ^V2{Ox) 

^^9i + i52)r(Ox)t, 
^^vtix) = T{Ox){vt{Ox) + dtMOxY , 
R : ■"■ui(x) = vi{Rx), = ~V2{Rx), 

'^vt{x) = vtiRx), 
T : ^v,ix) = t{Tx){v,{Tx) + d^)T{Tx)\ 

^vt{x) = ~T{Tx){vt{Tx) + dt)r{Tx)\ (3) 

where 1?^, with i e {1,2}, are the displacements along 
primitive translation vectors which arc chosen to be 

fli = (|a, -^a) and 02 = (0,V3a), respectively. Here 
a is the lattice spacing. Further, O, i?, and T in eq.([3]) 
represent a 60 degrees spatial rotation around the cen- 
ter of a hexagon, a spatial rcfiection, and time reversal, 
which are given by 

Ox ^ 0(:Xi,X2,t) ^ {\xi- ^X2,^Xi + \x2,t), 
Rx ^ R{xi,X2,t) ^ {Xi,-X2,t), 

Tx = T{xi,X2,t) = {xi,X2,-t), (4) 

respectively. In expressing these symmetry transfor- 
mation properties, we have introduced the matrix t{x) 
which takes the form 

r{x)^( 0. . -^"P^-*^(^))V (5) 

^ ' y cxp[iip(x)) J ^ ' 

Finally, the Abelian "gauge" transformation 

h{x) = exp(ia{x)a3) (6) 

belongs to the unbroken U{l)s subgroup of SU{2)s and 
acts on the composite vector fields as 

^U^y = vl{x)~df,a{x), 

v^i^y = «±(x)exp(±2za(z)). (7) 
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FIG. 1: Bipartite non-Bravais honeycomb lattice consisting of 
two triangular Bravais sublattices. 




FIG. 2: Brillouin zone of the honeycomb lattice with corre- 
sponding hole pockets. 



Analytic calculations as well as Monte Carlo simula- 
tions in i-J-like models on the honeycomb lattice have 
revealed that at small doping holes occur in pockets cen- 
tered at lattice momenta fc" = — fc^ = (0, J^l^ ), and 
their copies in the periodic Brillouin zone |62l . |63| . The 
honeycomb lattice, illustrated in figure 1, is a bipartite 
non-Bravais lattice which consists of two triangular Bra- 
vais sublattices. The corresponding Brillouin zone and 
the corresponding hole pockets are shown in figure 2. 
The single-hole dispersion relation for the t-J model on 
the honeycomb lattice is illustrated in figure 3. 

The effective field theory is defined in the space-time 
continuum and the holes are described by Grassmann- 
valued fields tp[{x) carrying a "flavor" index / = 
that characterizes the corresponding hole pocket. The 
index s = ± denotes spin parallel (+) or antiparallcl (— ) 
to the local staggered magnetization. As will be shown 




FIG. 3: Energy-momentum dispersion relation Eh{k)/t for 
a single hole in the t-J model on the honeycomb lattice for 
J/t = 2. 



in [65|, under the various symmetry operations the hole 
fields transform as 



SU{2)s : 


ipiix)' 


= exp{±ia{x))ipj^{x), 


U{1)q ■■ 




= exp{iuj)ipj_{x) , 


: 


^•t^iix) 


= exp{ik{ ai)ip^{x) , 


O : 


o^jlix) 


= zf exp{TMOx) ± i^)^p'^{Ox), 




O^^lix) 


= T cxpiTi^iOx) T if 


R : 


^^Plix) 


= i^i{Rx), M(^) =V'±(i?a;), 


T : 


^^lix) 


^cxp{Ti^{Tx))i,i\Tx), 






^ cxp{TiLp{Tx))il;'^\Tx), 






= — exp{ztiip{Tx))ijjj.{Tx), 




^^i\x) 


= -exp{±iip{Tx))^l{Tx). (8) 



Here U{1)q is the fermion number symmetry of the holes. 
Interestingly, in the effective continuum theory the loca- 
tion of holes in lattice momentum space manifests itself 
as a "charge" k{ under the displacement symmetry Di. 

Once the relevant low-energy degrees of freedom have 
been identified and the transformation rules of the cor- 
responding fields have been understood, the construc- 
tion of the effective action is uniquely determined. The 
low-energy effective action of magnons and holes is con- 
structed as a derivative expansion. At low energies, terms 
with a small number of derivatives dominate the dynam- 
ics. Since the holes are heavy nonrelativistic fermions, 
one time-derivative counts like two spatial derivatives. 
Here we limit ourselves to terms with at most one tem- 
poral or two spatial derivatives. One then constructs all 
terms consistent with the symmetries listed above. The 
effective action can be written as 

5[4^V'i,e1= / d'xdt ^£„^, (9) 

where denotes the number of fermion fields that the 
various terms contain. The leading terms in the pure 
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magnon sector take the form 



A = -4- [ die- die+ —dte- dte 



(10) 



The leading terms with two fermion fields (containing at 
most one temporal or two spatial derivatives) are given 
by 

/ = Q,/3 
s = + ,- 

+ A?/>f (iswj +cr/W^)V/, 
+ iK[{Di +isafD2)i}jl\vl + isa fvl)ij}[,^ 
- {v{ + isafv'2)^i\Di + isafD2)4'Ls] 

+ isafN2 (^Af ^wjw^"*^/ - ^4'' '"2^1 "^4) ■ 

(11) 

Note that all low-energy parameters that appear above 
take real values. It should be noted that vf{x) contains 
one spatial derivative, such that magnons and holes are 
indeed derivatively coupled. In eg. pTj) . M is the rest 
mass and M' is the kinetic mass of a hole, A is the leading 
and if is a subleading hole-one-magnon coupling, L, Ni 
and a-re hole-two-magnon couplings, and 

ff^{x)=d.v^ix)-d,v!ix) (12) 

is the field strength of the composite Abelian "gauge" 
field. The sign erf is + for f = a and — for f = (3. The 
covariant derivative in eq. (|lip takes the form 

D^^Piix) = d^^iix) ± ivl{x)i;i{x). (13) 

The leading terms with four fermion fields and without 
derivatives are given by 

^4 = E {^(cV.-^-lV'^^. + ^fVfV'^t^l) 

+ G2CV"V'f^f + G3C^CV'^l^^.} (14) 

with the real- valued 4- fermion coupling constants Gi , G2 , 
and G3. In principle, there are even more contact in- 
teractions among the fermions, such as 6- and 8-fermion 
couplings as well as 4-fermion couplings including deriva- 
tives. Since these terms play no role in the present work, 
we will not list them explicitly. 

Remarkably, the leading terms of the above Lagrangian 
have an accidental continuous 0(7) rotation symmetry 



FIG. 4: Homogeneous phase with constant staggered magne- 
tization. 

that acts as 

0(''Vi(a:) - exp(Ti7/2)V'i(0(7)^), 
'^^'^'^vi{x) — COS7 wi (0(7)2:) -|-sin7 W2(0(7)x), 
'^^'^'^V2{x) — — sin7 i;i(0(7)a:) +COS7 V2{0{'i)x), 
0(7)x = 0(7)(xi,X2,i) = 

(cos7 xi — sin7 X2, sin7 xi -\- COS7 X2, t). (15) 

This symmetry is not present in the underlying micro- 
scopic systems and is indeed explicitly broken by the 
higher-order terms in the effective action. 

III. HOMOGENEOUS VERSUS SPIRAL 
PHASES 

This section is devoted to the analysis of homogeneous 
and spiral configurations of the staggered magnetization, 
illustrated in figures 4 and 5, respectively. First, the 
energy of doped holes is calculated keeping the staggered 
magnetization field fixed. Then the parameters of the 
staggered magnetization field are varied in order to min- 
imize the total energy. 

A. Fermionic Contribution to the Energy 

In this subsection we compute the fermionic contribu- 
tion to the energy of a homogeneous or spiral configu- 
ration of the staggered magnetization. For the moment, 
we ignore the 4-fermion couplings. The considerations of 
this paper are valid only if the 4-fermion couplings are 
weak and can be treated in perturbation theory. Fur- 
thermore, we may neglect the vertices proportional to 
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was shown in (5j| , the most general configuration of this 
kind represents a spiral in the staggered magnetization. 
The Hamiltonian can then be diagonalized by going to 
momentum space. Since magnon exchange does not mix 
the flavors, the Hamiltonian can be considered separately 
for / = a and / = /?, but it still mixes spin s = + with 
s = — . The single-particle Hamiltonian for holes with 
spatial momentum p — (pi,P2) takes the form 



Hf{p) = 



A/ I 2A/' 

A(— ici + (JfC2) 



K{ici + cr/C2) 

2M' 



M 



(20) 



The hole-one-magnon vertex proportional to A mixes the 
spin s = + and s = — states and provides a potential 
mechanism to stabilize a spiral phase. The diagonaliza- 
tion of the above Hamiltonian yields 



FIG. 5: Spiral phase with helical structure in the staggered 
■magnetization. 



Eiip) 




(21) 



K, L, A^i, and which involve two spatial derivatives 
and are thus of higher order than the hole-one-magnon 
vertex proportional to A. The fermion Hamiltonian re- 
sulting from the leading terms of the Euclidean action is 
given by 



H 



J E [ 



1 



2M' 



The covariant derivative takes the form 



A:*£(a;) 



(16) 



(17) 



Here '^{'^ (x) and {x) are creation and annihilation op- 
erators (not Grassmann numbers) for fcrmions of flavor 
f = a,f3 and spin s = (parallel or antiparallel to 

the local staggered magnetization), which obey canoni- 
cal anticommutation relations. As before, aa = 1 and 
(T/3 = —1. The above Hamiltonian is invariant against 
time- independent U{l)s gauge transformations 



*£(x)' = exp(±ia(a;))^'£(a;), 

vfix)' = vf{x)-dMx), 
vf{xy = vf (x) exp{±2ia{x)) . 



(18) 



Here we consider holes propagating in the background 
of a configuration with 



(19) 



where cf and are real-valued constants. In other words, 
we focus on configurations of the staggered magnetiza- 
tion in which (after an appropriate gauge transformation) 
the fermions experience a constant composite vector field 
Di(x)', which leads to a homogeneous fermion density. As 



where |c| = c\ + c^- Interestingly, the above equation 
is independent of the flavor index /. We will keep the 
flavor index to indicate that there are two flavors in our 
calculations. Since the energy depends only on |c|, unlike 
in the square lattice case, potential spiral configurations 
do not prefer any particular spatial direction. This is due 
to the 0(7) spatial rotation symmetry discussed in the 
previous section. However, one should keep in mind that 
0(7) is an accidental symmetry of just the leading terms 
in the effective action, which is broken explicitly by the 
higher-order terms. Hence, when the higher-order terms 
are included, one expects the spiral to align with a lattice 
direction. Mixing via the A vertex lowers the energy e[_ 

and raises the energy E:!^_. It should be noted that, in this 
case, the index ± no longer refers to the spin orientation. 
Indeed, the eigenvectors corresponding to £"£ are linear 
combinations of both spins. The minimum of the energy 
is located at p = for which 



EiiO) = M + 



,3^2 



MI 

2Af' 



±A|c| 



(22) 



Since cf does not affect the magnon contribution to the 
energy density, we fix it by minimizing e[ (0) which im- 
C2 = 0. The energies of ea. (pT|) then reduce 



plies c- 
to 



Ei{p)=M 



2M 



-±A\c\. 



(23) 



Consequently, the filled hole pockets P^. are circles de- 
termined by 



2M' 



T. 



±1 



(24) 



f . . . 

where Ti^ is the kinetic energy of a hole in the pocket 
at the Fermi surface. The area of an occupied hole 
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pocket determines the fermion density as 



/ 



1 



(27r)2 Jpf 



1 

2^^ 



-M'T: 



Irpf 



(25) 



The kinetic energy density of a filled pocket is given by 



1 



(2^)2 J PL 



n 



1 



2M' 47r 



-M'T: 



The total density of fermions of all flavors is 



1 



and the total energy density of the holes is 



(26) 



M'{T^ + + + ), (27) 



(28) 



Here eo is the energy density of the system at half-filling. 
For 2'Kps > M'A^ the energy is minimized for Ci = and 
the configuration is thus homogeneous. The total energy 
density in the four-pocket case is then given by 



£4 = eo + Mn + 



4M' 



(34) 



For 2ttPs < M'A^, on the other hand, the energy is not 
bounded from below. In this case, |c| seems to grow with- 
out bound. However, according to eq. (|32p this would lead 



to Tjf < which is physically meaningless. What really 
happens is that two pockets get completely emptied and 
we are naturally led to the two-pocket case. Before turn- 
ing to that case, for completeness we first discuss the 
three-pocket case. 



C. Three Populated Hole Pockets 



with 



e£ = {M ± A\c\)ni + ti. 



(29) 



The filling of the various hole pockets is controlled by the 
parameters t£ which must be varied in order to minimize 
the energy while keeping the total density of holes fixed. 
We thus introduce 



S = £h- pn, 



(30) 



where ^ is a Lagrange multiplier that fixes the density, 
and we demand 



dS 



1 



j = —M'{M±A\c\+Tl-p) 



dTl 27r 



0. 



(31) 



B. Four Populated Hole Pockets 

We will now populate the various hole pockets with 
fermions. First, we keep the configuration of the stag- 
gered magnetization fixed and we vary the t£ in order 
to minimize the energy of the fermions. Then we also 
vary the parameters Ci of the staggered magnetization 
field in order to minimize the total energy. One must 
distinguish various cases, depending on how many hole 
pockets are populated with fermions. In this subsection, 
we consider the case of populating all four hole pockets 
(i.e. with both flavors f = a, (3 and with both energy 
indices ±). In this case, ea. (PT|) implies 



fi = M + 



T. 



f 



2M' ' ^ =^ 2Af ' 
The total energy density then takes the form 

e = eo + e,„ + eh 

= €o + 2pM' + (-+ + ^-+4 + ^- 

= eo + 2ps\c\' + Mn+ — --M'A^\c\\ (33) 

4M TT 



F7TA|c| 



(32) 



We now populate only three pockets with holes: the 
two pockets with the lower energies and as well 
as the pocket with the higher energy E". Of course, 
alternatively one could also flU the /3+-pocket. We now 
obtain 



eh = e"+e^+el 



(35) 



such that cq. f3T|) yields 



2Trn A , , 



27r7T, 4A , 



3M' 

rj--iCt _ 



27rn 2A, , 
3M^ + T"^"- 



(36) 



The total energy density then takes the form 

e = eo + e,„ + e,, = eo + 2p,|cp + e^ + e^ + e^ 

= eo + 2p4c\^+(M--\c\]n+ — 



2 

3^' 



(37) 



For Sirps > M'A^ the energy density is bounded from 
below and its minimum is located at 

, , TT An 

1^1 = 4 3^^^-W 

The resulting energy density in the three-pocket case 
takes the form 



es = eo + Mn 



3M' 



1 M'A2 
1 



n^. (39) 



•,3ttps-M'A^^ 

It is energetically less favorable than the homogeneous 
phase because 63 > €4 for 2n ps > M' A^ . For 2-Kps < 
M'A^ one obtains T" < which is unphysical. In fact, 
the Q;-)--pocket is then completely emptied and we are 
again led to investigating the two-pocket case. 
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D. Two Populated Hole Pockets 

We now populate only two pockets with holes. These 
are necessarily the pockets with the lower energies EZ 

and . In this case we have 

n = n^+n'! = ^M'(T^+T^), e,, = e^+/, (40) 
and thus cq. (|3ip now implies 

,-M+^-A\cl T, = T^^^. (41) 
The total energy density then takes the form 

e = eo + era + Cft = eo + 2ps|cp + el + e'^ 



eo + 2p,,|c|2 + (M- A|c|)n + 



TTn 



2M' 



(42) 



The energy density is bounded from below and has its 
minimum at 



4ps 

The value at the minimum is given by 



£2 = eo + }^^n + 



A2 



2M' 8p, 



(43) 



(44) 



The two-pocket spiral phase is less stable than the ho- 
mogeneous phase if £2 > £4, which is the case for 27rps > 
M'A^. As we have seen, both the three- and the four- 
pocket calculation become meaningless for 27rps < M'A^, 
because the kinetic energies Tjf then become negative 
which is unphysical. The two-pocket calculation, on the 
other hand, continues to make sense for 27rps < Af'A^. 



E. One Populated Hole Pocket 

Finally, let us populate only one hole pocket, say the 
states with energy E". Of course, alternatively one could 
also occupy the /?_-pocket. One now obtains 

T_"^|?. (45) 

The total energy density then takes the form 

9 

TTT) 

= eo + 2p,,|cp + (M-A|c|)n+— , (46) 



which is minimized for 



I I ^ 



(47) 



and the corresponding energy density takes the form 



ei = Co + Mn + 



(48) 



The one-pocket spiral is always energetically less favor- 
able than the two-pocket spiral. 



IV. INCLUSION OF 4-FERMION COUPLINGS 
IN PERTURBATION THEORY 

In this section the 4-fermion contact interactions are 
incorporated in perturbation theory. Depending on the 
microscopic system in question, the 4-fermion couplings 
may or may not be small. If they are large, the result 
of the perturbative calculation should not be trusted. In 
that case, one could still perform a variational calcula- 
tion. In this work we limit ourselves to first order per- 
turbation theory. We will distinguish four cases: the ho- 
mogeneous phase, the three-pocket spiral, the two-pocket 
spiral, and the one-pocket spiral. Finally, depending on 
the values of the low-energy parameters, we determine 
which phase is energetically favorable. 



A. Four-Pocket Case 

Let us first consider the homogeneous phase. The per- 
turbation of the Hamiltonian due to the leading 4-fcrmion 
contact terms is given by 



AH 



(49) 



It should be noted that ^'{^(x) and ^'{(a;) again are 
fermion creation and annihilation operators (and not 
Grassmann numbers). In the homogeneous phase the 
fermion density is equally distributed among the two spin 
orientations and the two flavors such that 



(50) 

The brackets denote expectation values in the unper- 
turbed state. Since the fermions are uncorrelated, for 
/ ^ /' or s ^ s' one has 



/t,T//\T;-'''t\Tc'''\ — /\T//t,T;/ 



(51) 



Taking the 4-fermion contact terms into account in first 
order perturbation theory, the total energy density of 
eq. ([M]) receives an additional contribution and now reads 



£4 = eo + 



4A-/' 



^{Gi + G2 + G^)n\ (52) 



Three-Pocket Case 



For a spiral aligned along the 1-direction (ci > 0, C2 = 
0) with cf = the eigenvectors of the single-particle 
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Hamiltonian of ea. pO) ) corresponding to the energy eigen- 
values Ej_ (p) are given by 



1 

V2' 



-i,(*^-*^:).(53) 



Inserting this expression in cq. (P5|) aUows us to evaluate 
the expectation value (AH) in the unperturbed states 
determined before. In the three-pocket case the states 
with energies E"{p), E^{p), as well as E"{p) (or alter- 
natively E^{p)), and with p inside the respective hole 
pocket are occupied and one arrives at 



1 



M'A 



I k2 



2 i-Kp, - M'k^ ) 3 ' 

1 M'A2 



4 Bttp, - i\/'A2 ) 3 • 
As a result, the energy density of eg. ([55)1 turns into 



(54) 



63 = eo + Afn + 
47rp, 



3A/' 
Af'A2 1 



1 

1 - - 



M'A2 



(37rp^ - Af'A2)2 32 
(4Gi + 3G2 + 3G3)Af 'A^ 



37rp,, -M'A2^ 
8(Gi + G2 + G3)^Ps 



(55) 



Two-Pocket Case 



In this case only the states with energy E"{p) and 

E^ {p) with p inside the respective hole pocket P[ are 
occupied and hence 



(56) 

As a result the energy density of eq. (P^ turns into 

= .o + Mn+ - n^ + i(G2 + G3).^. (57) 



D. One-Pocket Case 

In the one-pocket case only the states with energy 
E"{p) (or alternatively with E^{p)) and with p inside 
the corresponding hole pocket are occupied so that one 
has 

(58) 

In this case, the 4-ferniion terms do not contribute to the 
energy density which thus maintains the form of ea. (|48)) . 
i.e. 

.,^.o + Mn+(^-^)n^. (59) 




HomogGnGous phi 



homogeneous phase 



FIG. 6: The compressibilities Ki as functions of M' h? /2-Kps 
determine the stability ranges of the various phases. A ho- 
mogeneous phase, a spiral, or an inhomogeneous phase are 
energetically favorable, for large, intermediate, and small val- 
ues of ps, respectively. 



E. Stability Ranges of Various Phases 

Let us summarize the results of the previous subsec- 
tions. The energy densities of the various phases take the 
form 



1 2 

ei = eo -I- Mn -I- -K^n 



(60) 



According to eas. ((59)) . (|57)) . ([55]), and ([52)) . the compress- 
ibilities Hi are given by 



K2 = 



K3 



A^ 

4p/ 

" ^'+\g 



27r 



M' 
27r 



4p,, 



G3 



1 - 



1 



M'A 



/ a2 



Aiips - Af 'A2 1 r 



(37rp,, - A/'A2)2 16 
(4Gi + 3G2 + 3G3)A/'A2 

TT 1 



8(Gi + G2 + G3)7rp, 



K4 



2A/' 



r(Gi 



'3; 
G2 



G3). 



(61) 



The compressibilities functions of M' fs? jl-Kps are 

shown in figure 6. For large values of ps, spiral phases 
cost a large amount of magnetic energy and the homoge- 
neous phase is more stable. To be more precise, in this 
regime one has K4 < K3 < K2 < Notice that ki is al- 
ways larger than K2 for any value of p^. As p^ decreases 
and reaches the value 



A/'A2 
27r 



(Ar)2A2G' 

4^^2 



(62) 
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at leading order in the 4-fermi couplings one finds K2 = 
K3 = K4. For smaller values of Ps, the two-pocket spiral is 
energetically favored until K2 becomes negative and the 
system becomes unstable against the formation of spatial 
inhomogeneities of a yet undetermined type. 

It should be pointed out again that these results apply 
only if the 4-fcrmion contact interactions arc weak. Even 
if the 4-fcrmion couplings arc indeed small, the results 
presented in this work do not necessarily reveal the true 
nature of the ground state. Due to the variational nature 
of the calculation, one cannot exclude that the phases 
that we found may still be unstable in certain parameter 
regions. It is instructive to compare the results presented 
here with the results obtained in the square lattice case 
[s^ l . Qualitatively the stability ranges of various phases 
are the same for both lattice geometries except that the 
one-pocket spiral is never energetically favored on the 
honeycomb lattice while it is favorable in a small param- 
eter regime on the square lattice. 



we have seen that, for sufficiently large values of ps, 
the homogeneous phase is energetically favored. With 
decreasing ps, a two-pocket spiral becomes energetically 
more favorable. On the other hand, in contrast to 
the square lattice case, the one-pocket spiral is never 
favored. For small values of p^ the two-pocket spiral be- 
comes unstable against the formation of inhomogeneities 
of a yet undetermined type. In [63j the low-energy 
parameters p^, c, and M' have been determined in terms 
of the parameters t and J of the underlying t-J model. 
It will be interesting to also determine the strength of 
the hole-one-magnon vertex A in order to decide which 
phase is realized in this model. Further applications of 
the effective theory, including the one-magnon exchange 
potential and the resulting two-hole bound states, are 
currently under investigation. 
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